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Abstract 

New study confirms conclusions made in [6]; according to it, there is a 
disturbed region expended along the CME propagation direction in front 
of a coronal mass ejection whose velocity u is lower than the critical uc 
relative to the surrounding coronal plasma. The time difference brightness 
(plasma density) in the disturbed region smoothly decreases to larger 
distances in front of CME. A shock wave forms at u higher than uc in the 
front part of the disturbed region manifested as a discontinuity in radial 
distributions of the difference brightness. 

1 Introduction 

A coronal mass ejection (CME) structure in white light is often characterized by 
the following well-known features: a bright frontal structure (FS) that covers the 
region of decreased plasma density (cavity) that may includes a bright interior 
(core). However, besides the said features, another extended disturbed region 
defined by [I] can exist immediately in front of a CME. The aim of our study 
is to investigate changes in the disturbed region form, when a CME velocity 
increases, and possibilities for formation of a shock wave in this case. 

2 Method of analysis 

In the analysis, corona images obtained with the SOHO/LASCO C2 and C3 
[I] were represented as the difference brightness AP = P(t) — P(to), where 
P(to) is the undisturbed brightness at to before the event considered, Pit) is 
the disturbed brightness at any instant t > to- Calibrated LASCO images 
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Figure 1: a, b - the slow CME1 on 5 May 1997; c, d - the fast CME2 on 20 
September 1997; top panels (a, c) display images of the difference brightness 
in the form of isolines; the bottom panels (b, d) present distributions of the 
difference brightness depending on the distance r, measured from the frontal 
structure center (point O) along two different sections "a" and "b", whose di- 
rections are shown on the top panels. 

were used with the total brightness P(t) expressed in units of the mean solar 
brightness (P ms6 ). 

Images of the difference brightness were employed to investigate the CME 
dynamics and disturbed region. For the purpose we used presentations in the 
form of isolines and sections both along the solar radius at fixed position angles 
PA and non-radial sections at various instants t. On all the images, the position 
angle PA was calculated counterclockwise from the Sun's north pole. 

3 Data analysis 

First we consider two CMEs (CME1 and CME2) whose velocities V differ greatly 
at R = 4-5 R Q (R© is the solar radius). Two top panels of Figure 1 show 
the typical difference brightness form (in isolines) for these two CMEs at the 
instants, when their frontal structures FS appear in the C2 field of view. Fig- 
ure la presents the slow CME1 (5 May 1997, V fa 400 km s" 1 ), Figure lc 
the fast CME2 (20 September 1997, V « 800 km s" 1 ). The values of ve- 
locity V correspond to the linear fit of the height-time measurements for the 
fastest frontal parts of the CMEs and were taken from the CME catalogue 
( http: //cdaw.gsfc. nasa.gov/CMEJist/[ ). 
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Figures la,b show that on the images of both CMEs the frontal structure 
FS can approximately be presented by a part of a circle with its center at O 
(dots on the figures) . The main direction of the CME propagation that roughly 
coincides with its symmetry axis is indicated by a heavy dashed line "b". It 
passes along the streamer belt or streamer chains [3J[5]; i.e., it is in the region 
of the quasistationary slow solar wind (SW). 

In order to find the left boundary of the disturbed region (from the CME 
side), by analogy with [TO] we determine the FS width h as a width at a half- 
height of the difference brightness AP(r) distribution constructed from the FS 
center. For CME1, the frontal structure in the direction of section "a" (Fig- 
ure la) is least distorted by the disturbed region effect and has a minimum 
width h « 0.3 R0 (a curve with light circles in Figure lb). Let us take the 
FS right-hand boundary as the left-hand boundary of the disturbed region, as 
shown in Figure lb. Its position is indicated by a vertical dashed line. In Fig- 
ure lb, a curve with solid circles shows the AP(r) distribution along the section 
"b" in the direction of the CME propagation. The disturbed region is marked 
by a horizontal line with an arrow-head and labeled respectively. 

The comparison between CME1 (slow) and CME2 (fast) yields two principal 
distinctions: 

1 . The isolines that correspond to the minimum difference brightness of the 
slow CME1 are extended along the direction of its propagation, while 
those of the fast one are close in form to a circle. 

2. The difference brightness AP(r) distribution along the direction of the 
CME propagation continuously decreases up to the most specific front 
part of the disturbed region for the slow CME, whereas in the front part 
of the fast CME disturbed region a discontinuity appears in the AP(r) 
distribution on a scale dp « 0.2-0.3 R Q (shaded in Figure Id). 

The AP(r) distributions constructed along the direction of the CME prop- 
agation are presented for the set of eight CMEs in Figure 2. CME velocities are 
different and increase from bottom to top in Figure 2, thus the slowest CME 
with V ~ 230 km s _1 is shown on the bottom panel, and the fastest one with 
V ~ 2000 km s _1 on the top panel. 

Slanted hatching in Figure 2 a-h presents the disturbed region and half of 
the frontal structure (labeled as FS). As the CME velocity increases from min- 
imum (V w 230 km s _1 ) to critical one Vc ~ 750-800 km s , width of the 
hatched region tends to decrease. Discontinuity in AP(R) distributions at the 
front boundary of the disturbed region (shaded parts in Figure 2 a-e) is ob- 
served when CME velocities V are higher than the critical velocity Vc- At the 
same time at V < Vc such discontinuity is absent, and the difference brightness 
distribution smoothly decreases with increasing distance until it becomes indis- 
tinguishable. Notice that CMEs with velocities V close to the critical velocity 
VC (Figure 2 d-g) propagate near the plane of the sky; thus, their measured 
velocity is close to the real one (Figure 2 shows coordinates of places of the 
CME initiation on the solar disk, taken from [2 and Solar- Geophysical Data 
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Figure 2: The AP(R) distributions for 8 CMEs along the direction of their 
propagation. 
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(http://sgd.ngdc.noaa.gov). Appearance of such a discontinuity at the front 
boundary of the disturbed region at V > Vc was first described in [6] . 

Obviously, in processes of the "CME - undisturbed coronal plasma" inter- 
action a crucial role should be played not by a V value, but a value of CME 
velocity relative to the surrounding SW stream u — V — Vsw- Since CME veloc- 
ities were determined in the direction of their propagation, we took a velocity of 
the slow SW flowing for the most part in the region of the coronal streamer belt 
and streamer chains, along which the majority of CMEs move, as the velocity 
Vsw of the undisturbed solar wind [5] . 

Figure 3 presents values of the relative velocity u measured for twenty four 
different CMEs at different distances. For Vsw{R) we employed dependence 
derived by [12] of the slow SW velocity on the distance R in the streamer belt. 
This dependence is shown by a dash-dot line in Figure 3. 

In Figure 3, solid marks correspond to the CMEs having a discontinuity in 
the difference brightness distributions in front of the disturbed region. CME 
velocities V were determined from the discontinuity motion. Open marks in 
Figure 3 indicate the CME without discontinuity. In this case we took a velocity 
from the CME catalogue (http://cdaw.gsfc.nasa.gov/CMEIist/). Figure 3 shows 
that the cases with the discontinuity observed are in the high-velocity region, 
and the cases without discontinuity (the disturbed region smoothly decreased 
with distance is observed there) are for the most part in the low- velocity region. 
Hence we can assume that the discontinuity forms, when the relative CME 
velocity u exceed some critical uc value. A critical velocity value may depend 
on a distance R. 

Compare the obtained uc value with the typical velocity of disturbance 
propagation in the magnetized corona plasma that is roughly equal to the fast- 
mode MHD velocity Vms m the plasma. In order to estimate Vms at these 
distances we may use Va assuming that Vms ~ Va- Dashed line in Figure 3 
indicates the Va (R) dependence obtained by [9] . 

Obviously in Figure 3 the Alfven velocity passes approximately between 
clusters of points, which apply to the CMEs with discontinuity and without it. 
Hence uc ~ Va, i.e., the desired critical velocity is roughly equal to the typical 
velocity of disturbance propagation in the magnetized plasma. 

Hence we have a situation the classical gas dynamics refers to as "transonic 
transition" and formation of a shock wave. It was predicted theoretically, but 
it is first observed experimentally in the magnetized plasma. 

4 On possibility for resolution of a shock front 
width 

Here we will briefly discuss the problem of a possibility for resolution of a shock 
front width in the corona. The discontinuity is observed in distributions of 
the difference brightness AP(R) that results from free-electron scattering and 
is averaged along the line of sight in the optically thin corona. Since we do 
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Figure 3: CME velocities u = V — Vsw relative to the surrounding SW depend- 
ing on a distance from the solar center for the CME without discontinuity (open 
marks) and the CME with discontinuity (solid marks). The dash-dot curve in- 
dicates the velocity Vsw 01 the quasistationary slow SW in the streamer belt 
from [T2]. The dash curve shows the Alfven velocity in the streamer belt from 
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not know exactly the matter-density distribution along the line of sight, the 
observable scale Sf in AP(R) distributions may differ from a real scale 5n of 
the plasma density discontinuity. As a result of the averaging the observable 
discontinuity in the difference brightness profile can have larger scale than the 
real discontinuity in the density profile has. 

In order to estimate an effect of such averaging, Sf/Sn ratios were found 
in the context of a simple geometrical shock- front model in [6j. In the model 
considered, the shock-wave front was represented as a spherical shell with an 
outer radius Rf] the center of the shell was in the plane of the sky at Rc from 
the solar center (these parameters were specified according to the CME form). 

The brightness distribution AP(R), induced by free-electron scattering within 
the shell in the range from the shell center to its front edge, was calculated. At 
the given distance R, the brightness value is defined by the integral along the 
line of sight: 



where i(R,9) is the brightness induced by the one-electron scattering, N(r) - 
density. The i(R,0) function depends on a impact distance R and an angle 6 
relative to the plane of the sky. In the spherical shell, the density was supposed 
to change only depending on the distance r from the sphere center. In each case 
we chose a density profile N(r) such that a model brightness profile AP(R) 
obtained from N(r) by integration of Equation (1) showed the best correlation 
with the experimental profile of the difference brightness AP(R). Then a scale 
S n of the plasma density discontinuity was found from the density profile N(r) 
obtained. 

Calculations show that the AP(R) profile broadening does not exceed 20% 
in comparison with the density profile N(r). Thus, scale of a brightness pro- 
file discontinuity is a good approximation for determining scale of a density 
discontinuity in a shock wave front. 

It is worth noting that, at R > 10 Rq, a new discontinuity (with thickness 
Sp <C Sp) is observed to form in the anterior part of the shock front for events 
with V > Vc considered above. Within the experimental error, thickness S F ss 
0.1-0.2 Rq does not vary with distance and is determined by the LASCO C3 
instrument spatial resolution. Transfer from the shock front with thickness Sf to 
the discontinuity with thickness S F <C Sf can be interpreted as transition from 
the collisional to collisionless shock wave [7] . Similar discontinuities in brightness 
profiles associated with collisionless shock waves were registered ahead fast halo- 
type CMEs (V > 1500 km s" 1 ) at distances > 10 R in p]. 

5 Conclusions 

It has been shown that in front of a coronal mass ejection having a velocity u 
lower than the critical uq relative to the surrounding coronal plasma there is 
a disturbed region expended along a direction of the CME propagation. The 




(1) 



7 



time difference brightness AP in the disturbed region smoothly decreases up to 
larger distances in front of the CME. Given u > uc, a discontinuity forms in 
distribution of difference brightness or plasma density in the disturbed region 
front part. Since the uc value is close to the local fast-mode MHD velocity, 
which in corona approximately equal to the Alfven one, the formation of such 
a discontinuity when uc is exceeded may be identified with the formation of a 
shock wave. 
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